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Propagation  of  first-order  electromagnetic  discontinuities 

in  an  isotropic  medium. 

* 

by 

R.  Venkataraman 
and 

R.  S.  Rivlin 

Center  for  the  Application  of  Mathematics 
Lehigh  University 
Bethlehem,  Pa. 

Abstract 

The  propagation  of  a  first-order  electromagnetic  disconti¬ 
nuity  is  discussed.  Expressions  are  obtained  for  the  possible 
velocities  of  propagation  as  functions  of  the  field  strengths 
ahead  of  the  surface  of  discontinuity.  Expressions  are  also 
obtained  for  the  growth  in  the  magnitude  of  the  discontinuity 
as  the  wave  progresses. 


-2- 


1 .  Introduction. 

In  this  paper  we  shall  discuss  the  propagation  of  electro¬ 
magnetic  waves  in  which  the  fields  are  continuous  with  respect 
to  space  and  time,  but  where  spatial  and  time  deriatives  may  be 
discontinuous  on  some  surface.  Such  discontinuities  are  called 
first-order  discontinuities.  Until  recently  studies  of  the 
propagation  of  discontinuities  of  this  type  have  been  confined 
to  fluid  dynamics  (see,  for  example,  Courant  and  Friedrichs^). 

More  recently,  using  a  technique  based  on  the  compatibility 
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conditions  at  a  moving  surface  of  discontinuity,  Thomas  has 
studied  the  growth  of  elastic  waves  of  this  type  in  linear  ma¬ 
terials.  Varley  and  Cumberbatchm  have  considered  the  more 
general  problem  of  a  quasilinear,  first  order,  system  of  hyper¬ 
bolic  equations  using  a  variant  of  the  method  for  a  linear 
system  described  in  Courant  and  Hilbert^.  In  this  paper  we 
apply  the  procedure  of  Varley  and  Cumberbatch  to  the  study  of 
the  propagation  of  first-order  electromagnetic  discontinuities 
in  a  non-linear  centrosymmetric  isotropic  material,  for  which 
the  electric  displacement  field  depends  only  on  the  electric 
field  and  the  magnetic  induction  field  depends  only  on  the  mag¬ 
netic  intensity  field. 

We  obtain  first  the  secular  equation  for  the  propagation 
of  a  first-order  electromagnetic  discontinuity  surface.  We 
find  that  there  are,  in  general,  two  possible  propagation 
velocities  in  a  specified  forward  direction  and  correspondingly 
two  possible  polarization  directions.  Both  the  velocities  and 
directions  of  polarization  depend  on  the  magnetic  induction  field  and 
electric  displacement  field  immediately  ahead  of  the  disconti¬ 
nuity. 

In  §4,  we  obtain  an  expression  for  the  growth  of  the 
magnitude  of  the  discontinuity  as  a  surface  of  discontinuity 
propagates.  In  §  §  5  and  6  this  result  is  specialized  to  two 
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cases  in  which  the  fields  ahead  of  the  surface  of  discontinuity 
are  uniform  and  constant.  In  one  of  these  the  magnetic  induction 
field  is  taken  normal  to  the  surface  of  discontinuity  and  in  the 
other,  the  electric  displacement  field  is  taken  to  be  normal  to 
this  surface. 


2. 
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Propagating  Surfaces. 

In  this  section,  we  consider  the  propagation  of  an  electro¬ 
magnetic  wave  for  which  the  magnetic  induction  field  B  and  ; 
elastic  displacement  field  D  are  continuous,  while  their  first 
derivatives,  both  with  respect  to  space  and  time,  may  be  dis¬ 
continuous  across  some  surface.  This  singular  surface  is  called 
a  first-order  singular  surface. 

Let  x  be  the  vector  position  of  a  generic  point  on  this 
singular  surface  with  respect  to  a  fixed  origin.  Let 


<(>(x,t)  =  0  , 


(2.1) 


for  fixed  t,  be  the  equation  of  the  singular  surface  at  time 
t.  We  define  a  by  : 

a  -  <p  (x,  t)  ,  :  (2.2) 

and  assume  that  <f>  is  a  continuous  function  of  x  and  of  t, 
which  has  continuous  first  derivatives  with  respect  to  space 
and  time.  We  also  assume  that  8<{i/3t  ^  0  in  some  neighborhood 
of  the  propagating  singular  surface,  i.e.  a  =  0. 

In  view  of  the  definition  (2.2)  of  a,  we  see  that 
a  =  constant  describes  a  propagating  surface,  which  is  the 
singular  surface  if  the  constant  is  zero.  Also,  the  variables 
x,a  define  a  point  on  such  a  surface.  Let  n(x,a)  denote  the 
unit  normal  at  the  point  x  to  the  surface  a  =  constant,  drawn 
in  the  direction  of  propagation.  Let  V(>0)  denote  the  speed 
of  propagation  of  the  surface.  From  the  definition  of  the 
speed  V  of  a  moving  surface  at  a  point  as  the  speed  with  which 
the  surface  moves  normal  to  itself  at  that  point,  we  have 


(2.3) 


n  • 


dx 

dt 
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where 


n  =  V<t>/(V<t>,-V<p)  x/2  . 


(2.4) 


Now,  from  (2.2)*  we  have  for  the  moving  surface  a  =  constant, 


dx 


!♦  '  31  +  H  -  0  ' 


(2.5) 


Multiplying  (2.5)  throughout  by  n  and  using  (2.3)  and  (2.4), 
we  obtain 


n  V<j> 

V  "  a*/ at 


(2.6) 


Now,  any  function  f,  say,  of  x  and  t  may  also  be  regarded 
as  a  function  7,  say,  of  x  and  a,  where  a  is  defined  by  (2.2). 

We  then  have 

f(x,t)  =  7(x,a)  .  (2.7) 


Whence 


3f  37  3a_  =  57 
9 1  5a  3 1  3a  3t  ' 

We  define  the  operator  V  as  the  spatial  gradient,  hold¬ 
ing  a  constant.  It  follows,  from  (2.7),  that 

vf  =  vT  =  vf  +  H  v<j>  . 

~  ~  ~  9a  ~ 


(2.9) 


With  C 2 . 6 )  and  (2.8),  this  yields 


Vf  =  vf  -  5.  —  ii=vf-  -  — 

~  ~  V  3a  3 1  ~  V  3t 


Thus,  we  may  take  the  operator  v  as 


w  =  n  +  2.  L_ 
l  l  V  3t 


(2.10) 


(2.11) 


Here  we  may  point  out  that  the  operator  V  is  a  surface 

operator  on  the  surface  <p  =  0  since  by  the  choice  of  the  new 

set  of  coordinates,  the  variables  involved  in  the  differenti- 

4  - 

ation  correspond  to  a  displacement  on  the  surface  . 
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3.  The  Secular  Equation. 

In  the  absence  of  free  charges.  Maxwell's  equations  may 
be  written  as 

VxE  +  3B/3t  =  0  , 

(3.1) 

VxH  -  3D/3t  =  0  , 

where  E  and  H  are  the  electric  and  magnetic  intensity  fields 
respectively. 

We  shall  consider  that  the  medium  in  which  the  wave  is 
propagating  is  isotropic,  cent rosymme trie  and  non- dissipative 
and  that  E  is  determined  only  by  D  and  H  is  determined  only 
by  B.  It  follows  that  the  constitutive  equations  take  the 
forms 

E  =  e(D-D)D  and  H  =  y(B-B)B  ,  (3.2) 

where  e  and  y  are  scalar  functions  of  the  indicated  arguments. 
Introducing  (3.2)  into  (3.1),  we  obtain 

Vex  D  +  eVxD  +  |§-  =  e'V(D.D)xD  +  eVxD  +  f|-  =  0 

and  (3.3) 

3D  3D 

VyxB  +  yVxB  -  =  y'v(B-B)xB  +  uVxB  -  =  0  , 

where  e'  =  3e/3(D*D)  and  y'  =  3y/3(B*B)  . 

We  bear  in  mind  the  known  vector  identity 


V ( D • D)  =  2  [(D-V)D  +  D* (VxD) ] 


(3.4) 


i 
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Substituting  from  (3.4)  and  (2.11)  in  (3.3)^,  we  obtain 
SD  3D  SB 

.  e2x  It  +  2£'CD-3f)CnxD)  i-  V  - 

=  2e'V[(D-V)DxD  +  {Dx(VxD)}*D]  +  eVfxD!.  .  (3.5) 

Similarly,  from  (3.3)2  we  obtain 

SB  SB  SD 

ynx  —  +  '2y' (B*-g~)  (nxB)  +  V  ^ 

=  2y’V[ (B-V)BxB  +  (Bx(VxB) }xB]  +  yVVxB j  .  (3.6) 

We  now  introduce  the  assumption  that  B  and  D  are  con¬ 
tinuous  at  the  moving  surface  <f>  =  0 .  It  follows  that  their 
tangential  derivatives  VB,  VD,  VxB  and  VxD  are  also  con¬ 
tinuous  across  the  surface. 

We  denote  by  b  and  d  the  jumps  in  SB/St  and  SD/St  as  we 
cross  the  surface  <f>  *  0  in  the  direction  of  n;:  thus 

b  =  [SB/St]  ,  d  =  [SD/St]  .  j  (3.7) 

With  (3.7)  we  obtain,  from  (3.5)  and  (3.6),  on  =  0 

-  I 

enxd  +  2e  '  (D-d)  (nxD)  -  Vb  .  -  0 

ynxb  +  2y ' (B *b) (nxB)  +  Vd  =0 


and 


(3.8) 


From  (3.8),  we  obtain 


n*b  =  0  and  n*d  =  0  .  (3.9) 

Thus,  the  jumps  b  and  d  are  tangential  to  the  surface  <j>  =  0 . 

Now,  eliminating  d  from  (3.8),  we  obtain  on  <J>  =  0 

yenx(nxb)  +  2y ' e (B • b) (nx (nxB) } 

+  2e ' (nxD) [yD- (nxb)  +  2y ' (B-b) {D- (nxB) } ]  +  V2b  =  0  .(3.10) 

We  introduce  the  notation 

B  =  B .  +  B  ,  D  =  D  +  D  ,  (3.11) 

~~t~n,~~t~n’  v 

where  B^  and  are  tangential  to  the  surface  <j>  =  0  and  Bn 

and  D  are  normal  to  it.  We  also  write 
~  n 

J  =  nxpt  .  (3.12) 

Using  (3.11),  (3.12)  and  (3.9)1  in  (3.10),  we  obtain 
(V2  -ye  )  b  -  2y»e(Bt.b)Bt 

- 2e 1 J [y (b • J)  +  2y * (B .  -b) (B . • J)  ]  =  0  .  (3.13) 

Forming  the  dot  product  of  (3.13)  with  B  ,  we  obtain 
[(V2-ye)  -  2p»e(Bt.Bt)  -  4y ' e ' (Bt • J) 2 ] (Bt*b) 

=  2ye  * (b- J) (B  -J) 


(3.14) 
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Forming  the  dot  product  of  (3.13)  with  J,  we  obtain 
[-2y*e(J-Bt)  -  4y  *e  '  (B^  J)  (J.  J)  ]  (B^b)  ! 

=  [(ye-V2)  +  2ye ' (J- J) ] (b-J)  .  (3.15) 

Eliminating  b  from  (3.14)  and  (3.15)  and  using  the  notation 

x  =  V2  -  ye  ,  (3.16) 

we  obtain 

x2  -  x  [2y'{e(Bt-Bt)  +  2e*(Bt-J)2}+  2ye’(J-J)] 

+  4yy'ee'{(J-J)(Bt-Bt)  -  (J-Bt)2}  =  0  .  (3.17) 

From  (3.16)  and  (3.17),  we  obtain  the  two  solutions  for  V2: 

V2  =  ye  +  y'[e(Bt.Bt)  +  2e’(J-Bt)2]  +  pe';(J*J). 

±{[y,{e(Bt‘Bt)  +  2e*(J-Bt)2}  +  ye’(J-J)]2 

-  4yy * ee ' [ (J * J) (B  *B  )  -  (J-B.)2]1/2  .  (3.18) 

^  <v  ~  b  ~  t  -V  -V  L 

V2  is  real  and  positive  if  and  only  if 

[y'{e(Bt.Bt)  +  2e  '  (J  •  Bt)  2  }  +  ye  •  (J .  J)  2  ]  2  ;• 

5  4yv'ee'[(J*J)(Bt*Bt)  -  (J-Bt)2]  (3.19) 

and  ye  >  0.  In  this  case  we  see  that  for  given  n,  B  and 
Dt>  there  are  two  possible  real  speeds  of  propagation  in  the 
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forward  direction.  Propagation  in  the  opposite  direction  with 
equal  speeds  is  also  possible. 

Provided  that  V2  is  a  simple  root  of  the  equation  (3.17), 
i.e.  the  greater  than  sign  in  (3.19)  is  valid,  we  can  determine 
the  direction  of  polarization  of  b  from  equation  (3.10).  Then 
the  direction  of  polarization  of  d  can  be  determined  from 
equation  (3.8)2* 
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4 .  Growth  Equation.  : 

In  this  section  we  shall  derive  a  differential  equation 
governing  the  change  in  magnitude  of  the  discontinuity.  In 
order  to  do  this  it  is  more  convenient  to  refer  all  equations 
to  a  rectangular  cartesian  reference  system  and  to  use 
indicial  notation.  With  this  notation,  we  obtain,  by  differ¬ 
entiating  equations  (3.3)  with  respect  to  t, 


3  3D£ 

ei;jk{Pk£  ft<5xj>  +  Q 


3D  3D0  3  2B .  : 

_H  _A}  +  _ L  =  o 

nk£  3t  3x .  a.2  * 

J  O  L 


(4.1) 


~  3B  9B  *  9B  0 

p  3  s  £  ^  ,  qi  n.  £ 

eiik  k£  9t^9x.^  nk£  3t  9x. 

J  J 


9  2D. 


-  k  0  , 


where 


Pvo  “  +  2e 1 DVD0  , 


+  2y  1  B,B Q  , 


=  3Pk£/3Dn  =  2e'Vk£  +  4e”DnDRD£  +  >2e'Dk6£n 


(4.2) 


+  2e’Dn6 


£  kn  » 


S.vo  =  3RkJl/3Bn  =  2U’B„6Vo  +  4p"Br>BvB0  +  -2y»Bv« 


n  k£ 


n  k  £ 


k  £n 


*  2«'Vkn-  • 


and  e ^ denotes  the  alternating  symbol. 
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We  now  change  the  independent  variables  from  x^t  to  xi,a 
In  indicial  notation,  the  relation  (2.10)  may  be  expressed  as 


3  -  r_i_)  .  11  !_ 

377  "  v  .at  ’ 

J  J  u 


(4.3) 


whence 


3  (  3  i  =  _3_r3_)  =  r  JL't  1_  -  -L  li¬ 
st  ^  3x  .  3x  .  ^  3 1^  ldx.J  at  V  9t2 
J  J  J 


(4.4) 


Introducing  (4.3)  and  (4.4)  into  (4.1),  we  obtain 


32D0  3  2B .  3  3D*, 

eijkPk£nj  Ft7”  "  V  ^T~  =  Veijk{PkJlCn7)a  3t“ 


9Dn  (C!V 

V-3x.-,a  v  at  r 


n .  3D, 


^nkJl  3t 


(4.5) 


32B0  32D.  9  3B^ 

eijkRk*nj  +  V  iT7"  =  Veijk{Rk£(axJ)o  Tt 

3B  /  3B „  n.  3 BA 

+  SnkJl  3t"(  (377^o  ‘  V1  3t")}  ’ 


Eliminating  V  32D^/3t2  from  (4.5),  we  obtain 


32B  3D 

(ais+V25is^=  eijk{v2pMC5lT)at3tt) 


+  v2Q 


nki  3 1 


n  l  c^1)  _  1± 

t  V  l3x.J  V  3 1  / 


(4.6) 
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VeS,mnPk)lRnsn3  ')a  ^t 

•'  m 


/  9B  9B 

'Ve£mnPkJlSpnsnj  V^d)a  ^  ' 


n  9B  9B 
m  _ s_  p 

V  9t  9t  J 


L>  • 


where 


a.  =  e.  .,e.  P,  .R  n  .n 
is  13k  &mn  kH  ns  j  m 


(4.7) 


Using  square  brackets  to  denote  the  jump  in  a  quantity  across 
the  surface  of  discontinuity,  we  have 


9Dn  9D* 


9D_  9D, 


9D  9Dn 
nw  £■, 


+  csr'icw4)  -  .  (4.8) 


where  (9Dn/9t)Q  denotes  the  value  of  9Dn/9t  immediately  in 
front  of  the  wave.  Also,,  we  note  that  P^,  ^nk£’  Rk&  and 
Snk£  are  continuous  across  the  surface  of  discontinuity. 

Using  these  results,  we  obtain  from  (4.6),  with  the 
notation  b^  =  [9B^/9t]  ,  d^  =  [9D^/9t], 


32Bq  9dp 

tais+v2{is)t^]  -  aijk{V2Pk,(Hj)a 


9D0  n.  / 

+  y2Qnkil^'93c7)  dn  '  (dnd£ 
3  a  v 

9D  9D  ,v 

*  <?A  dn  +  W1)  ^ 


-  Ven  P,  „R  n  .  (tt— =-) 
Jlmn  kS,  ns  3  ^  9x  J 

J  TO 
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3B  n 

VeAmnPkASpnsnj  ^3xm-*a  bp  '  V~  *-bsbp 


3B  3B 

+  CH2-)o  bs  +  <ar>0  V}} 


We  now  introduce  the  notation 


bs  *  ors  ’ 


where  r  is  a  unit  vector  in  the  direction  of  b  .  With 
s  ^ 

(4.10)  and  (4.2),  equation  (3.8)2  becomes 
VdS,  =  'RnseAmnnmbs  “  "RnseAmnnmrsa 

Using  this  result  to  substitute  for  Vd^  in  (4.9),  we 
obtain, 


32B 


^is'^is^TfT4  =  Li  +  MiCT  *  Nia  ’ 


where 


L. 

x 


_  r  f  do  •.  r  3a  ■.  i 

-Ve.  e0  P,  0R  t  (n  (-5— — )  +  n-  (3-5—)  ■>  > 

ljk  Amn  kA  ns  s  m  3x^  a  1  a 


Mi  '  "y2eijkeAmnPkA('3x^a  Rnsnmrs^ 

3DA 

+e  i  j  kQnk  ARtpns  rp { “  ens  tV  ( 3x7^  a 

3Dp  3Dn 

+ens tnj  ('3t-)0  +  eAstnj  ^Ft-^  } 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


(4.13) 


'Vei  jkeAmnPkS,nj  {Rns  ('Sx") 

•'  »  tn 


3B  n  3B  n  8Bo 

*W^a  rP  ‘  r1  rs ■rV^»  ’ 


P  eijkenste3lab^nkJl^tp^bcnjnansrpr( 


teijketanPMSpnsnjnmrsrp  ' 


In  our  present  notation  (3.8),  may  be  written 


^i  _  eijkPUnjdt 


With  (4.11),  we  obtain 


(4.14) 


(a.  +V26 .  )b  =  0  , 

v  IS  IS  s  ’ 


(4.15) 


where  a.  is  given  by  (4.7).  We  define  a  unit  vector  H •  by 
1  s  1 


^(ai^V^is)  -  0 


(4.16) 


Multiplying  (4.12)  throughout  by  % ■  and  using  (4.16),  we 
obtain 


L  +  Ma  +  No2  =  0 


(4.17) 


where  L,  M  and  N  are  defined  by 


(L,M,N)  =  £i(Li,Mi,Ni) 


(4.18) 
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Introducing  (4.7)  and  (4.10)  into  (4.15)  and  differ¬ 
entiating  with  respect  to  n^,  we  obtain 


(eijkelmnPkiRns<6pjV6pmV  +  2V5^6is)rs 


+  ^ei  jkekmnPkkRnsnjnin+^  ^is^3n  ~  ^ 

P 


(4.19) 


Multiplying  throughout  by  and  bearing  in  mind  (4.16)  and 
(4.7),  we  obtain 


ijkeS,mn^kS.^ns  ^pjnm+(^pmnj^irs 


2VlrrVi 

P 


(4.20) 


Since  a^g ,  defined  by  (4.7),  is  homogeneous  of  second 
degree  in  n^  and  V  is  determined  by  the  discriminant  of 
equation  (4.15),  it  follows  that  V  is  homogeneous  of  first 
degree  in  n^.  Therefore, 


V  =  n  9V/3n 
P  P 


(4.21) 


Now,  suppose  a  point  x^  depends  on  t  in  such  a  manner  that  it 
always  lies  on  the  moving  surface  of  discontinuity.  Then 
its  velocity  is  dx^/dt.  The  velocity  of  the  surface  normal 
to  itself  is  V,  which  is  therefore  given  by 


V  =  n  dx  /  dt 
P  P 


(4.22) 


From  (4.21)  and  (4.22),  we  obtain 


n  3V/3n  =  n  dx/dt 

p  p  P  P 


(4.23) 
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A  possible  choice  for  the  dependence  of  x^  on,  t  which  satisfies 
the  condition  that  x^  shall  always  be  on  the  surface  of  dis¬ 
continuity  is  that  for  which 


dx  /dt  =  3V/9n 
P  P 


(4.24) 


We  shall  make  this  choice*  and  calculate  the  growth  in  the 
magnitude  of  the  singularity  as  we  move  along  the  locus  of 

V 

From  (4.13)^,  (4.18),  (4.20)  and  (4.24),  we  obtain 


L  -  2V*  lihHr)  Vi  ■  2V2  &  ^iri 

p  pa 


(4.25) 


We  introduce  (4.25)  into  equation  (4.17)  and  make  the 
transformation 


a  =  1/u  . 


(4.26) 


We  obtain 


a!  -  Mu  +  N  , 


(4.27) 


where 


M  =  M/2V2£iri  ,  '  N  =  N/2V25-iri 


(4.28) 


This  equation  has  the  general  solution 


u  -  i  -  e'Udt 
c 


{c  +  /n  e"/Mdtdt  }  , 


(4.29) 


*In  this  case  the  locus  of  x  is ,  in  fact,  a  bicharacteristic 

p 

of  the  secular  equation  |  sl .  g V2 ^ ^ s  I  =0.  ; 
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where  C  is  a  constant  of  integration. 

We  see  that  if  the  fields  B  and  D  are  known  before  passage 
of  the  surface  of  discontinuity,  the  velocity  of  such  a  surface, 
propagating  in  a  given  direction,  can  be  determined  from  (3.18). 

and  Rj^  defined  by  (4.2)  can  be  determined  and  hence  a^s  can 
be  determined.  r^  can  then  be  calculated  from  (4.15)  and  (4.10), 
and  l ^  from  (4.16).  Then,  M  and  N  can  be  determined  from  (4.28), 
(4.18)  and  (4.13). 
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5 .  Special  case:  magnetic  induction  field  ahead  of  wave  in 
direction  of  propagation. 

We  now  assume  that  the  fields  ahead  of  the  moving  surface 
of  discontinuity  are  uniform.  This  implies  that  the  speeds  of 
propagation  V  given  by  (3. 18)  are  constant .  It  then  follows, 
from  (4.15) ,  that  b  has  a  constant  orientation  (i.e.  r  is  a 
constant)  and,  from  (4.16),  that  the  unit  vector  £  is  constant. 
We  see,  from  (4.13),  that  in  this  case  =0  and  is  con¬ 

stant.  Then,  from  (4.18),  M  =  0  and  N  is  constant.  We  then  ob¬ 
tain  from  (4.27)  or  (4.29) 


u  ■  S'  '  Nt  *  5~  >  l5-1) 

o 

where  oQ  is  the  value  of  o  when  t  =  0.  The  jump  magnitude 
a  increases  or  decreases  accordingly  as  N  is  negative  or 
positive. 

In  the  particular  case  when  the  direction  of  propagation 
is  in  the  1- direction  and  the  applied  magnetic  induction  field 
into  which  the  wave  propagates  is  also  in  the  1-direction,  we 
have 


n-  =  6  . ,  and  B.  =  B<5 .  .. 

i  ll  i  ll 


(5.2) 


Introducing  (5.1)  into  (3.18),  we  have  since  B  =  0, 


V  =  ye  or  Vz  =  ye  +  2ye'(Dz+Dz)  . 


(5.3) 


Corresponding  to  the  first  of  these  results,  we  obtain 


r  =  (D2+D2)-‘/2  (0,d  .D.) 


(5.4) 
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and  corresponding  to  the  second,  we  have 


X  =  (D2+D2)"l/2(0, 


D  ,  D  ) 

3  ’  2 


(5.5) 


The  directions  of  polarization  of  the  b  vectors  associated 
with  the  two  waves  are  perpendicular  to  each  other. 

Again,  for  the  wave  corresponding  to  the  first  of  the 
solutions  (5.3),  we  find  that  N  =  0,  so  that  the  discontinuity 
propagates  unchanged. 

We  now  consider  the  wave  corresponding  to  the  second  of 
the  solutions  (5.3).  Introducing  (5.2)  into  (4.2)2,  we  obtain 


R 


ns 


=  y  6 


ns 


2y 1 B2 


^ni5  si 


(5.6) 


From  (4.2)^,  (4.7),  (5.2)  and  (5.6)  we  see  that  a^s  is  sym¬ 
metric  for  interchange  of  i  and  s.  Hence,  with  (5.5), 

*i  -  ri  -  (^+D23)-l/!C-D3«i2.D2S.3)  .  (5.7) 


From  (4.18)  and  (4.13),,  we  obtain,  using  (5.2), 


N  =  -V  e  ■  .e  .  en  ,Q  ,  „  R .  R.  r  r  l  - 
11k  nit  Jlib^nkS,  tp  be  p  c  1 


+  e-  1  e„  P,  nS  r  r  £. 
11k  Jtin  kfl,  pns  s  p  1 


(5.8) 


Introducing  the  expressions  (5.7)  for  r^  and  SL^  into  (5.8) 
and  using  (5.6),  (4.2)  and  (5.2),  we  obtain,  with  (5.3)2  and 
(4 . 2  8)  2 

jr  =  _  y2 (3e 1  +2e"D2)  D 
(ye+2ye ' D2) 5/ 2 


(5.9) 
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where  D  =  (Dj+D2)1/2  is  the  magnitude  of  the  component  of  the 
electric  displacement  field  normal  to  the  direction  of  propa¬ 
gation.  (In  carrying  out  the  calculation,  we  find  that  the 
second  term  on  the  right  hand  side  of  (5.8)  is  zero.)  In  thb 
case  when  D  =  0,  the  two  waves  with  velocities  given  by  (5.3) 
have  the  same  speeds  and  can  be  propagated  with  arbitrary 
direction  of  polarization  in  the  23-plane.  ■ 
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6 .  Special  case:  electric  displacement  field  ahead  of  wave 
in  direction  of  propagation. 

We  again  assume  that  the  fields  ahead  of  the  moving  surface 
of  discontinuity  are  uniform,  so  that  (5.1)  still  applies.  We 
again  take  the  1-direction  as  the  direction  of  propagation,  but 
now  assume  that  the  electric  displacement  field  is  in  the  1- 
direction,  so  that 


n.  =  6-  and  D.  =  D6-  .  (6.1) 

in  l  n  v 

It  follows  from  (3.12)  and  (6.1)  that 

J  =  0  .  (6.2) 

Introducing  this  result  into  (3.18)  we  obtain 

V2  =  ye  or  V2  =  ye  +.2y»e  B  *Bt  .  (6.3) 

Corresponding  to  the  first  of  these  results,  we  have 

r  =  (B2+B2)  _1/2(o,-B3,B2)  (6.4) 

and  corresponding  to  the  second  of  these  results,  we  have 

r  =  (B2+B2r1/2(0,B2,B3)  .  (6.5) 

From  (6.4)  and  (6.5)  it  follows  that  the  directions  of 
polarization  of  the  electric  vectors  for  the  two  waves  are 
perpendicular. 

For  the  wave  corresponding  to  the  first  of  the  solutions 
(6.3),  R-  =  0,  so  that,  from  (5.1),  the  discontinuity  propagates 
unchanged. 
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We  now  consider  the  wave  corresponding  to  the  second  o£ 
the  solutions  (6.3).  As  in  §5,  we  again  obtain  the  expression 
(5.8)  for  N.  Introducing  (6.1)  into  (4.2)^  wje  obtain 


2e '  D2  6. 

k  i 


(6.6) 


Proceeding  in  a  manner  analogous  to  that  used  in  the  previous 
section,  we  find  that 


B(3y  1  +2ii"B2 ) 
ye  +  2y'eB2 


(6.7) 


2  2  l/2  - 

where  B  =  (B2+B3)  >  is  the  magnitude  of  the  component  of  the 
magnetic  induction  field  normal  to  the  direction  of  propaga¬ 
tion.  (In  carrying  out  the  calculation,  we  find  that  the  first 
term  on  the  right  hand  side  of  (5.8)  is  zero.) 
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